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Nuclei close to the neutron drip line are described employing an interaction model which is based 
on the low-momentum interaction Viowk- This effective two-body interaction which is determined to 
reproduce the nucleon-nucleon (NN) scattering data at energies below the pion thresh-hold is sup- 
plemented by a density-dependent contact interaction fitted to reproduce the saturation properties 
of infinite nuclear matter within the Hartree-Fock approach. It is demonstrated that corresponding 
calculations for closed shell-nuclei using this interaction model reproduce the bulk properties of 
these nuclei, independent whether the wave functions are expanded in terms of harmonic oscillator 
waves or in a basis of plane waves discretized in a spherical box of appropriate size. This plane 
wave basis, however, is more appropriate to describe weakly bound nuclei and the transition from 
discrete nuclei to homogeneous matter which is supposed to occur e.g. in the crust of neutron stars. 
Properties of exotic nuclei are studied within a Hartree-Fock plus BCS approximation. 

PACS numbers; 21.60.Jz, 21.30.Fe, 21.65.-f, 26.60.Gj 

I. INTRODUCTION 

The new generation of radioactive beam facilities, e.g. the future GSI facility FAIR in Germany, the Rare Isotope 
accelerator planned in the United States of America or SPIRAL2 at GANIL/France, facilitate the nuclear structure 
studies away from the line of /3-stability, especially for the neutron-rich nuclei. The study of these nuclei is of high 
interest since they represent important steps in the nuclear reaction chains for the formation of elements. The weak 
binding energies of these nuclei and the corresponding large spatial extensions of the orbits for the valence nucleons 
lead to very interesting features like e.g. a neutron halo made of several neutrons outside the core. 

These nuclei can furthermore be understood as a first step in the transition from isolated nuclei to infinite homo- 
geneous matter, which should occur in the outer crust of a neutron star. This transition is a very good example 
for a quantum liquid with a phase transition from the droplet to the homogeneous phase. This transition is a chal- 
lenge for theoretical nuclear structure physics as it incorporates the transition from isolated nuclei via a crystal-like 
structure of quasi-nuclei embedded in a sea of neutrons to the phase of homogeneous baryon matter. These struc- 
tures have been described by means of Thomas-Fermi calculations or mean-field calculations in a Wigner-Seitz cell 
using sirnple phenomenological models for the nucleon-nucleon (NN) interaction like the density-dependent Skyrme 
forces[ll,i,i,i,[i,@] or a relativistic mean- field approach Q. 

All these investigations are based on phenomenological interactions, which have been adjusted to describe the 
saturation properties of symmetric nuclear matter and the structure of stable nuclei with large binding energies located 
in the valley of /^-equilibrium. These models and interactions have been developed to provide a simple description 
of the mean field in terms of local single-particle densities. The predictive power of such simple phenomenological 
interactions, however, may be rather limited as these models have been constructed to fit the data of finite nuclei. 

An alternative to this phenomenological approach is the use of so-called realistic interactions like the charge- 
dependent Bonn interaction 8] or one of the Argonne interaction models[9]. Such realistic interactions are based e.g. 
on the meson exchange model and adjusted to describe the experimental data of nucleon-nucleon (NN) scattering. 
These interactions contain rather strong short range and tensor components, which make it inevitable to employ non- 
perturbative approximation schemes for the solution of the many-body problem for the nuclear hamiltonian based on 
such interactions [13] ■ Nuclear structure calculations using such interactions have become feasible for very light nuclei 
and provide a very good description if an appropriate three-nucleon force has been added[llj. In the next future such 
sophisticated calculations will not be possible for heavier nuclei or the nuclear structures in the crust of neutron stars 
mentioned above. 

A possible way out of this problem is to consider an interaction model, which separates the low-momentum (below 
a cut-off A) and high- momentum components of a realistic NN interaction by means of renormalization techniques [l^, 
[H, 0, |T^. If the cutoff A is chosen around A = 2 fm~^ the resulting low-momentum interaction Viowk still describes 
the NN scattering data up to the pion threshold and turns out to be independent of the underlying realistic interaction 
V . Since the high-momentum components, which correspond to the short-distance behavior, of V have been removed, 
the resulting Viowk does not produce short-range correlations in a significant way, but can be treated within the 
Hartree-Fock approximation 15]. 

Employing Viowk in a calculation of nuclear matter, however, one obtains a binding energy per nucleon increasing 



2 



with density in a monotonic way, such that no saturation point is obtained[l5|, |16|- In order to compensate this 
deficiency we will supplement Viowk by a simple density-dependent contact interaction adjusted to reproduce the 
empirical saturation property of symmetric nuclear matter. This density-dependent contact interaction, which is 
introduced to obtain the empirical saturation point of nuclear matter, can also be interpreted as a parameterization 
of the three-nucleon force, which, as discussed above, seems to be necessary for a microscopic description of nuclear 
structure data based on realistic NN interactions. Therefore we had the hope to obtain a good description of finite 
nuclei using this interaction model. We will demonstrate below that this hope has been fulfilled to a good extent. 

Because of the non-locality of Viowk nuclear structure calculations of finite nuclei using this interaction are typically 
not performed in coordinate space but employ a set of appropriate single-particle wave functions as a set of basis 
states, typically those of a harmonic oscillator potential. The choice of the harmonic oscillator (HO) basis seems 
to be rather plausible if one wants to describe the structure deeply bound, double magic nuclei like ^7?e, ^^O, and 

Co [TtI. ITsj. This choice may be questionable if nuclei close to the proton or neutron drip line are to be considered. 
The HO basis may not be appropriate to describe the tail of the single-particle wave functions for the weakly bound 
valence nucleons. The main reason, why this oscillator basis is so popular, must be attributed to the fact that the HO 
basis provides a simple transformation from relative coordinates in coordinate or momentum space to a coordinate 
system located in a laboratory frame. Such a transformation is required to allow for structure calculation of finite 
nuclei using realistic interaction as defined e.g. by Viowk- 

The HO basis is certainly not appropriate to describe the quasi-nuclear structures in the pasta phase of neutron 
stars mentioned above. For those studies a basis of plane wave (PW) states confined to a finite box with a size large 
to the nuclear dimension may be more suitable and could be considered as an appropriate basis for a Wigner-Seitz 
cell calculation!^, 0]. One of the main aims of this work is to establish the techniques, which are needed to allow for 
nuclear structure calculations of finite nuclear systems using realistic interactions in such a PW basis. 

In this work we describe double magic nuclei and weakly bound neutron-rich nuclei in the framework of Hartree- 
Fock (HF) plus BCS calculation using a Viowk potential supplemented with a contact interaction. We are going to 
investigate nuclei up to ^°^P6 and compare the results of calculations using HO and PW basis states. 

The plan of this paper is as follows. The procedure to renormalize the low-momentum interaction and obtain the 
matrix elements of Viowk is reviewed in Sec. [Hi In this section we also describe our interaction model supplementing 
Viowk by an appropriate contact interaction. We also work out how to solve the HF equation in a basis of harmonic 
oscillator and plane wave states. This section is completed by a short description of the BCS procedure. The results 
for infinite matter and finite nuclei are presented and discussed in Sec. IIIII Finally, Sec. IIVI contains a summary and 
the conclusions of our work. 

II. MODEL OF THE NN INTERACTION 

The idea behind the Viowk interaction model is to separate the low-momentum and high-momentum components 
of realistic interactions and restrict the nuclear structure calculation to the low-momentum components. This means 
that one tries to define a model space, which accounts for the low-momentum degrees of freedom and renormalizes the 
effective Hamiltonian for this low-momentum regime to account for the effects of the high-momentum parts, which are 
integrated out. A model space technique to disentangle these parts, the unitary-model-operator approach (UMOA) 
[T9| . can be used to calculate the matrix elements for Viowk- 

First, we define the projection operators P, which projects onto the low-momentum subspace, and Q, which projects 
onto the complement of this subspace, as in all well known model space techniques. Furthermore, these operators 
satisfy the usual relations like P + Q = 1, = P, — Q, and PQ = = QP. The idea of the unitary-model- 
operator approach is now to define a unitary transformation U so that the transformed Hamiltonian does not couple 
P and Q, which means 

QU-^HUP^O (1) 
has to be fulfilled. The effective interaction is defined in terms of this unitary transformation as 

Veff ^U-\ho+Vi2)U -ho, (2) 

with vi2 representing the bare A^A^-interaction and a starting Hamiltonian ho describing the one-body part of the 
two-body system. The unitary operator U can be written as 

[/ = (l+w-cjt)(l + cjcj^+w'fw)-i/2 (3) 

with an operator w fulfilling uj — QujP and lu^ = uj^^ = 0. In [l^ the operator w is obtained by first solving the 
two-body eigenvalue equation 
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and afterwards defining the matrix elements of lo using the eigenstates |$p) having the largest overlap with the 
P space. With the corresponding U then the effective interaction Vg// is calculated as described in [l^j US]- In 
this way one obtains an effective Hamiltonian H^ff — ho + Vef / , which contains this effective interaction Vef / . The 
eigenvalues, which are obtained by diagonalizing this effective Hamiltonian in the model-space (P-space), are identical 
to the diagonalization of the original Hamiltonian ho + V in the complete space. 

This model-space scheme can now be applied to the effective two-nucleon problem by considering for the basis states 
of the 2-nucleon system the states identified by the relative momentum, its modulus and the corresponding partial 
waves. For a given partial wave the states of the P-space are identified as those states with a relative momentum 
smaller than a cutoff A. Applying the technique described above leads to the effective interaction Viowk ■ Solving now 
the Schroedinger equation for Viowk within the model space, e.g. by solving the Lipmann Schwinger equation for NN 
scattering with the cutoff A, yields the same phase shift as solving the Lipmann Schwinger equation without cutoff 
for the original interaction vi2. If the original interaction vi2 is realistic in the sense, that it has been fitted to the 
experimental NN phase-shifts, these phase shifts, up to the cutoff A, will also be reproduced by Viowk- 

Furthermore, the resulting Viowk is found to be essentially model independent. This means that the potential is 
independent of the underlying realistic interaction vi2, if the cutoff A is chosen around A = 2 fm~^. Thus, one is 
able to extract a low-momentum potential Viowk, which in a model independent manner describes the low- momentum 
component of realistic NN interactions in a more or less unique way. Note, that, by construction, this Viowk interaction 
is defined in terms of matrix elements in a basis of NN states labelled by relative momentum. 

The fact that the high-momentum or short-range components of realistic NN interactions have been integrated out 
by means of the unitary transformation of Eq. ([3]), seems to lead to problems in the description of the saturation 
behavior of nuclear matter. Because of this lack of short-range correlation effects, which are modified in the medium, 
the emergence of a saturation point is prevented in calculations of symmetric nuclear matter (Tsl . IT^ . In order 
to achieve saturation in symmetric nuclear matter, three-body interaction terms or density-dependent two-nucleon 
interactions are needed. Therefore, the effective interaction Viowk is supplemented by a simple contact interaction, 
which we have chosen following the notation of the Skyrme interaction to be of the form 

AV = AVo + AVs, (5) 

with 

AVo = ^<o[(2 + xo)p^ - (2x0 + l){pl + pI)] (6) 

AV3 = ^i3P"[(2 + x^)p' - (2x3 + l){pl + pDI (7) 

where pp and pn refer to the local densities for protons and neutrons while the matter density is denoted as p = pp + pn- 
The parameters of the contact interaction are ^Oi 2^0, ^3, X3 and a. As described below we have chosen a fixed value of 
a = 0.5 and xq = 0.0 and fitted to and ts and X3 in such a way that HF calculations using Viowk plus the contact term 
of Eq. ([5]) yield the empirical saturation point for symmetric nuclear matter and reproduce the symmetry energy at 
saturation density. 

Within the HF approximation this interaction model leads to a single-particle Hamiltonian for protons and neutrons 
(j/ = n,p) of the form 

HhF,v — Hkin + Hiowk,u + ^Hct.p + HcoulSup, (8) 

where H^in is the kinetic part and AHct.v originates from the contact interaction of Eq. This part is given by 

AHot,. = to/2[(2 + xo)p - (1 + 2xo)p,] + i3/24[(2 + X3){2 + a)p^+" 

-(2x3 + l){2p>. + ap'--\pi + pI)}]. (9) 

Furthermore, for the charged particles, the protons, one has an additional contribution to the hamiltonian, the 
Coulomb contribution. It is given by 



and 



where the direct term is 



-^Coul — Ucoul,dir U(Joul^exch, 

(10) 



U -Aize'^l Idr'r'^Pp{r')/r for r' < r 
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and the exchange term is 

Ucoul,exch = ( - ) pl'^- (12) 



3,1/3 



Note that AHct,i and Hcoui are locaL They are defined in terms of the single-particle densities pi resulting from 
the eigenstates of Hhf/i, which implies that they have to be determined in a self-consistent way, as usual for HF 
calculations. 

As it has been mentioned above, the effective interaction Viowk is non-local and defined in terms of matrix elements. 
This implies that the HF calculation has to be performed in a Hilbert space using an appropriate basis |a>,|/3>,.... 
The HF Hamiltonian is then expressed in terms of the matrix elements between these basis states < a\HHF\l3 > and 
the HF single-particle states |5'„ > are defined in terms of the expansion coefficients in this basis 

|*n >= ^ \a >< a\^n >= ^c„ala > ■ (13) 

a a 

The part of the HF Hamiltonian originating from Viowk can be expressed in terms of two-body matrix elements by 

< a\Hiowk\/3 >= X! < ailViowklPS > p^s , (14) 

7,(5 

with pys representing the single-particle density matrix. 

If we restrict the HF variational procedure to a spherical description of nuclei, as we will do in the following, the 
HF Hamiltonian is diagonal in the angular momentum quantum numbers j,l,rn and the expansion is simplified to 

00 N 

< f|*„ >= *„ijm(r) = ^ Cnilj^iljmir) « ^ Cnilj^iljm{r) , (15) 
2=1 i=l 

where Cniij are expansion coefficients and ^iijm{r) =< r|a > are the wave functions of our orthogonal basis. The 
number of basis states N has to be chosen to guarantee that the results are not affected by this limitation. Note, 
that using this coordinate space representation of the basis states allows the calculation of the matrix elements of 
< a\HHF\P > for local terms in Hhf in a straight forward way. 

In this paper calculations are done using two different orthogonal basis systems. These basis systems are respectively 
that of a free particle and that of a spherical harmonic oscillator. In these cases the wave functions of our orthogonal 
basis can be separated in a radial part and an angular part, 

^ujm (r) = (r^iljm) = Ra {r)yij,^ (z9, ^) . (16) 

In both basis systems, yijm represent the spherical harmonics including the spin degrees of freedom by coupling the 
orbital angular momentum I with the spin to a single-particle angular momentum j. 

The radial part Ru is different for the two basis systems. The first basis system considered is that of a particle 
moving free in a spherical cavity with a radius i?, the plane wave (PW) basis [5[. The radial wave functions Ru are 
then given by the spherical Bessel functions, 

Ra{r) ^ N,iji{kar), (17) 

for the discrete momenta ku, which fulfill 

Ra{R) - Na]i{kaR) = => ji[kaR) = 0. (18) 
The normalization constant in Eq. (jl7p is given by 

{ i-K^ iovl = Q, 

{ Ji-i(kuR) 

It ensures that the basis functions are orthogonal and normalized within the box, 

d^r<i>*ij^lr)<^^'Vfm'{r) = 5u'5iv5jf5mm' ■ (20) 
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Next, the basis system of the harmonic oscillator (HO) is considered. The radial part of the basis functions takes 
the form ^ 



a3 T{i + l + i/2) 



(r/aY exp(-0.5 r^a^) L'+'/'(r7a2), 



(21) 



where is a Laguerre polynomial, F is a gamma function, and a = ^JJiJimuj) is the oscillator length. These 

basis functions ^ujm are orthogonal and normalized within a box with an infinite radius. However, in practice one 
can use a radius R >> a due to the exponential factor exp(— 0.5 r'^ /a^) in the radial basis function of Eq. ([2T|l . 

The matrix elements of Viowk or other realistic NN interactions, which are required in Eq. (|14p . are easily calcu- 
lated within the HO basis: In a first step one can calculate the oscillator matrix elements using the momentum 
representation of the relative basis. The second step, transforming these matrix elements from the relative system to 
the laboratory s^tem is also straight forward using the well known transformation brackets of the Talmi-Moshinsky 
transformation [il [H ^ . 

While the corresponding first step, the calculation of matrix elements of Viowk in the relative basis is trivial in the 
case of PW states, the transformation from relative to laboratory coordinates is much less convenient in the PW than 
in the HO basis. The evaluation and use of these vector brackets has been described e.g. in 24, 25] and will be used 
here. 

The Hartree-Fock approximation discussed so far is a good approximation for double closed shell nuclei. In order 
to consider nuclei with open shells as well, we are going to use a rather phenomcnological treatment of pairing 
correlations. The pairing is described by a BCS approach with a smooth cutoff. For the BCS potential V{r) we have 
taken a density dependent effective interaction of zero range as introduced by Bertsch and Esbensen[26l|. 



Vir) = Vo [l.O-v ipir)/po)% 



(22) 



where po = 0.16 fm~^ is the saturation density. Furthermore, the values Vq = 481.0, rj = 1.0, and j3 = 1.0 are taken 
for the parameters Vo, 77, and f3. These parameters have been determined by Garrido et aL[27|. 
The anomalous density has the form of 



nljm 



1) TT —^nljmirY 



where Unijm and Vnijm are variational parameters solved from the BCS equations. 



nljm 



1 



^nljrn 



nljm 



nljjn 



1 



nljm 



(23) 



(24) 



with Enijm — \J {^rdjrn — ^fY + ^liijm energy of the quasi-particle state and tnijm the single-particle energy 

determined in the HF equations. Furthermore, '^nijmi'r) is the HF orbital wave function of Eg . [TSl and Wnijm is the 
smooth cutoff, 

Wnljra = [1 .0 + exp((e„;j,„ - Ep - -Bcut)/^^)] (25) 

with Ecut = 5 MeV and Ae — i?cnt/10. Therefore, the local gap function can be written as 

A(r) = V{r)x{r), (26) 
where V{r) is the BCS potential of Eg. and x{t) is the anomalous density in Eq. (^5)) . 



III. RESULTS AND DISCUSSION 



A. Nuclear Matter 



All the results which we are going to discuss in this chapter employ a low- momentum interaction Viowk, which is 
based on the proton-neutron part of CD Bonn potentialQ]. Using a cutoff parameter A = 2 fm~^, these results do 
not significantly depend on this choice of the underlying realistic interaction. 

First, let us turn to the binding energy of symmetric nuclear matter, which are displayed in Fig. [1] The HF 
calculation using Viowk does not lead to a minimum in the energy versus density plot (see solid line in Fig.[T]), as we have 
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FIG. 1: (Color online) Binding energy per nucleon of symmetric nuclear matter of a Hartree-Fock calculation using Viowk (solid 
line) and of a calculation supplemented with a contact interaction (dashed line). Furthermore, comparison is made with a Dirac 
Brueckner Hartree Fock (DBHF) calculation [§3, [s^] (dashed-dotted line) and a Brueckner Hartree Fock calculation using the 
potential Argonne V18 plus a 3-nucleon force (BHF + 3BF) ^ (dotted line). 



already mentioned before [T^, [T^. The absence of saturation is one of the major problems of calculations for nuclear 
matter using Viowk- This problem cannot be cured by the inclusion of correlations beyond the HF approximation [T^. 
Using a conventional model for a realistic interaction, like e.g. the CD Bonn interaction or the Argonne V18, which 
are not reduced to its low-momentum components, one obtains a saturation point in calculations which account for 
correlations beyond the mean-field approximation. Therefore one may argue that the Viowk approach cannot reproduce 
the saturation of nuclear matter as it misses the quenching of short-range correlations in the nuclear medium, which 
is included in sophisticated calculations using one of the conventional models for a realistic NN interaction. 

Note, however, that calculations like Brueckner Hartree Fock (BHF) for these interactions lead to a saturation point, 
but are not able to reproduce the empirical data. Depending on the amount of tensor and short-range correlations 
produced by these interactions either the predicted density is too large or the calculated binding energy is too 
weak IC , [2g| . In order to obtain the experimental data for the saturation point one has to include a three- body 
force 29 or include relativistic effects e.g. within a Dirac Brueckner Hartree Fock (DBHF) approach [sol [Sll. [s^. l33j. 
The relativistic effects due to the change of the nucleon Dirac spinors in the medium must be considered within the 
framework of non-relativistic calculations by means of a density-dependent two-nucleon or an effective thrcc-nuclcon 
interaction. 

In Fig. m our calculation for the binding energy is compared to the ones of a DBHF calculation from Ref. [H, [s^ 
and the non relativistic BHF calculation from Ref. [l^ . The non relativistic BHF calculation is based on the Argonne 
Vis and includes a three-body force deduced from the meson-exchange current approach [29| . 

Using the Viowk approach, the absence of saturation can be cured by adding three-body forces [l3j . Therefore, we 
also supplemented the Viowk potential with the contact interaction of Eq. ([5]) to obtain the experimental saturation 
behavior with a binding energy per nucleon of E/A = —16.0 MeV at a saturation density po=0-16 fm~^. It is worth 
to notice that only the isoscalar part of the contact interaction can influence the binding energy in symmetric nuclear 
matter. 

In order to fix the isovector part of the contact interaction as well, we also inspect the symmetry energy of infinite 
matter as a function of density in Fig. O The symmetry energy of 21 MeV at saturation density obtained from the 
Viowk calculation is only about two third of the experimental value. Furthermore, the value of the symmetry energy 
is significantly below the one predicted by DBHF also at higher densities. This means that the symmetry energy, 
respectively the isospin dependence of the EoS, of Viowk without contact interaction is in principal too soft. 

So we have added a contact interaction term as defined in Eq. ([5]) choosing a value for a —0.5 and xq — 0.0. The 
three remaining parameters have then been fitted to reproduce the empirical saturation point of nuclear matter and 
the symmetry energy at saturation density. The results for these fitting parameters are listed in Table U and the 
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FIG. 2; (Color online) The symmetry energy as a function of the density. The symmetry energy of the Hartree-Fock calculations 
using both Viowk (solid line) and Viowk plus a contact interaction (dashed line) is plotted. For comparison the results of the 
DBHF approach [33] are plotted as well. 
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TABLE L Parameters io, xq, ta and X3 defining the contact interaction of Eq. ((SJ 

corresponding results for the binding energy of symmetric nuclear matter and the symmetry energy as a function of 
density are displayed in Figs. [T] and [2] as well. 

The first observation that becomes evident from Fig. [T] is that the Viowk plus contact interaction (Viowk+CT) and 
DBHF calculations yield results which are rather similar, while the result of the BHF approach predicts lower energies 
at high densities. The compressibility modulus of Viowk+GT at saturation density is K = 258 MeV, while those of 
the other microscopic calculations range from 200 MeV to 265 MeV. This means all the presented equations of state 
(EoSs) in Fig. [1] can be characterized as rather soft, at least at densities up to about three times saturation density. 
This prediction of a soft equation of state (EoS) is also supported from observables extracted from heavy ion reactions. 
For example, heavy ion data for transverse flow [3^ or from kaon production [ssj support the picture of a soft EoS in 
symmetric nuclear matter. 

In Fig. [2] we can see that the addition of the contact interaction enables us to bring the value of the symmetry 
energy close to the experimental value. The symmetry energy at saturation density is now 33 MeV. Moreover it is 
comparable to that of the microscopic DBHF approach up to densities of three times saturation density. Therefore 
the isospin dependence of the EoS for Viowk+CT calculation is rather stiff. 

B. Finite Nuclei 

In the first part of this section we are going to consider the bulk properties of nuclei with closed shells for protons 
and neutrons. The typical example for a heavy nucleus is the example of ^^^Pb. Results for HF calculations using 
the Viowk interaction, expanding the single-particle wave functions in a basis of harmonic oscillator (HO) states are 
presented Fig. [31 The left panel of this figure displays the calculated binding energy as a function of the oscillator 
length a (see Ea. (l2T]) ). The expansion of the single-particle wave functions has been restricted to oscillator states with 
N < 8. One finds that the resulting energy is rather sensitive to the choice of the oscillator parameter and leads to 
very attractive values at small values of a. This gain in binding energy is accompanied by a decrease of the calculated 
radius of the neutron distribution as can be seen from the right panel in this figure. The nucleus is collapsing to a 
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FIG. 3: (Color online) Results of Hartree Fock calculations for the energy per nucleon (left panel) and the radius of the neutron 
distribution (right panel) of ^"^Pb. The calculations use either the bare Viowk or the Viowk plus contact (CT) interaction. The 
single particle wave functions are expanded in an oscillator basis. Results are presented as a function of the oscillator length 
parameter a. 




1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 
Oscillator length [fm] 



FIG. 4: (Color online) Results of Hartree Fock calculations for the energy per nucleon of ^®0. The calculations use either the 
bare Viomk or the Viomk plus contact (CT) interaction. The single particle wave functions are expanded in an oscillator basis. 
Results are presented as a function of the oscillator length parameter a. 



tightly bound system with a large density. So these results reflect the properties of infinite matter if the bare Viowk 
interaction is employed. 

The situation is rather different if we consider light nuclei. As an example we focus our attention on ^^O and 
present in Fig.[4]the results for the energy per nucleon using various values for the oscillator parameter a. Using Viowk 
one obtains results which are rather insensitive on the choice of a over a large range of values. Only a basis with a 
larger than 2 fm seems not to be adequate for such calculation and leads to smaller values for the binding energy per 
nucleon. Also it is worth noting that the calculated energies around -7.7 MeV per nucleon are rather close to -7.98 
MeV observed in experiment. This rath er g ood description of light nuclei using Viowk in an appropriate HO basis has 
been observed before by Coraggio et a/. jlSl. ISGt- 
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FIG. 5: (Color online) Comparison of total density distributions calculated for ^^O using plane wave basis in a spherical box 
of radius ii = 15 fm (PW, solid lines) and an appropriate harmonic oscillator (HO, dashed lines) basis. Results of HF Yio-wk 
calculation without contact term (dark lines) and with contact term (light lines) are shown. 

How are these results for finite nuclei modified by the contact interaction, which we have added to describe the 
properties of infinite matter? This is visuahzed in Figs. [3] and |4] as well. In the case of ^°*P6 the contact interaction 
provides a stabilization of the results leading to values for the energy per nucleon and radius of the nuclcon distribution, 
which are in good agreement with the experimental data. 

The same is true as well for the case of ^^O. The results obtained for the energy per nucleon in HF calculation is 
very insensitive to the oscillator parameter considered and the calculated energy is rather close to the experimental 
value. 

The calculated energies for ^^O with and without the contact interaction are rather close to each other (see Fig. 2]), 
the resulting wave functions, however, are quite different. This can be seen from Fig. [5] showing the calculated nuclear 
density profiles. The bare Viowk interaction yields densities in the center of the nucleus, which are close to three times 
the saturation density of nuclear matter. Consequently, the calculated radii are very much below the empirical result. 
This is in line with the HF calculations of 36||. This also explains why oscillator parameters a larger than 2 fm, are 
not suitable to describe such a solution. The density profile derived from the Viowk + CT interaction model yields 
densities below the saturation density and consequently results for the radius of the charge distribution, which are in 
good agreement with experiment (see discussion below). 

Up to this point we have restricted the discussion to HF calculations performed within an oscillator basis. How do 
these results compare with the ones derived within the PW basis. Before we proceed to this comparison it should be 
mentioned that the calculations within the HO basis are much easier to perform than those within the PW basis. Not 
only that the techniques to calculate the two-body matrix elements (in particular the transformation from relative 
to laboratory coordinates) is more tediously, also the number of basis states which have to be considered is typically 
larger in the PW basis as compared to the HO basis. In order to obtain results which are stable from the numerical 
point of view, we have typically considered a reference box with a radius i? which is about 5 times as large as the 
radius of the nucleus considered. In the case of ^^O we have typically considered boxes with radii R — Vo fm, while 
a calculation of "^^^Pb seems to require a box as large as i? = 30 fm. This seems to be necessary in order to obtain 
a set of radial momenta {ku in Eq. pT]) ). which has a sufficient resolution. For small values of the reference box we 
also run into the problem that one is considering interactions between neighbored Wigner Seitz cells. For such large 
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Nucleus 




E/A [MeV] 






[fm] 




[fm] 




HO 


PW 


exp. 


HO 


PW 


HO 


PW 




-7.77 


-7.91 


-7.98 


2.68 


2.68 


2.66 


2.66 




-8.33 


-8.57 


-8.55 


3.41 


3.41 


3.37 


3.37 




-8.13 


-8.42 


-8.67 


3.46 


3.45 


3.68 


3.68 




-7.42 


-7.75 




3.60 


3.59 


4.11 


4.15 


208 


-7.25 


-7.76 


-7.87 


5.48 


5.45 


5.76 


5.70 



TABLE II: Results for the binding energies per nucleon and rms radii for different nuclei using the Viowk plus contact term 
interaction model using the harmonic oscillator (HO) or plane wave (PW) basis. The calculated energies have been corrected 
by subtracting the spurious energy of the center of mass motion. The experimental values are taken from Ref. [37| . 









16^ 












Neutron 






Proton 




Orbital 


HO 


PW 


exp. 


HO 


PW 


exp. 


Sl/2 


-37.024 


-37.162 


-47 


-33.478 


-33.601 


-44 ±7 


P3/2 


-19.730 


-20.006 


-21.839 


-16.358 


-16.632 


-18.451 


Pl/2 


-16.217 


-16.484 


-15.663 


-12.888 


-13.155 


-12.127 


^5/2 


-3.509 


-3.739 


-4.144 


-0.444 


-0.690 


-0.601 


lSl/2 


-1.288 


-1.566 


-3.273 


1.436 


0.839 


-0.106 


^3/2 


0.720 


0.339 


0.941 


3.562 


1.886 


4.399 



TABLE HI: Single particle energies for the orbital levels of O nucleus derived from HP calculations using the Viowk plus 
contact interaction. Results are displayed using the harmonic oscillator (HO) and plane wave (PW) basis. The experimental 
values are taken from Ref. [36l ]. 

boxes, however, we need a large number of such momenta to cover the relevant range of momenta. This number is 
typically of the order 15 or larger. The number of different oscillator states, which are required to obtain stable results 
is significantly smaller. 

Comparing the results displayed by dashed and dotted lines in Fig. [S] one find that the results for the density distri- 
butions, which are obtained on the PW basis (solid lines) are very close to the ones derived in the HF approximation 
using the HO basis. This is particularly true for the calculations including the contact interaction term. Note that 
the HF calculation using the PW basis leads to an energy per nucleon of -7.91 MeV, which is slightly more attractive 
than the corresponding number of -7.77 MeV obtained in the HO basis. This indicates that the larger number of 
variational parameters, which are optimized in the PW basis (see discussion above) leads to deeper minimum than 
the HF variational procedure in the HO basis. Both results are rather close to each other and to the experimental 
value (see results listed in table HI)) . 

This table exhibits very similar findings for the HF calculations of other nuclei with closed shells for protons and 
neutrons including examples from light (^^O) to heavy nuclei (^"^Pfo). It seems that our fitting procedure for the 
contact term in infinite matter yields very good results for the energies and radii of finite nuclei. 

Also the single-particle energies are very similar using the HO or PW basis. This can be seen from table lllll which 
lists the single-particle energies for the example of ^^O for the occupied Os and Op shells as well as for the IsOd shell. 
The single-particle energies evaluated in the PW basis are slightly more attractive than those derived within the HO 
approximation. The difference is smaller of the tightly bound Os states than for the more weakly bound states in the 
Op and IsOd shells. 

The shapes of corresponding single-particle wave function displayed in Fig. [H] may exhibit the reason for this 
behavior. The wave functions resulting from the HO and PW cannot be distinguished in this logarithmic plot for 
radii up to 7 fm. For larger distance from the center, however, the oscillator model yields a faster decrease reflecting 
the Gaussian asymptotic, which is typical for the HO model. 

The differences in the asymptotic behavior of the wave functions are larger, if we consider the shapes of the most 
weakly bound orbits for protons and neutrons in the case of the weakly bound ^°Ca, which are shown in Fig.[71 This 
is in line with our expectations that the PW basis is superior to the HO basis in particular for the description of 
weakly bound nucleons. 

Finally, we want to study the effect of the basis to be used on the predictions for the occurrence of the neutron 
drip line. For that purpose we perform HF plus BCS calculations, as described in section [Hi using the HO as well as 
the PW basis. As an example we present the results for Fermi energy of neutrons for the A^e isotopes ranging from 
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FIG. 6: The neutron wave functions for the state pi/2 states of ^''O calculated in the HF approximation using two different 
basis systems: the plane wave basis (solid lines) and the harmonic oscillator basis (dashed lines). 

A = 20 to v4 = 32 in Fig. [H One can see that the differences in the predictions are rather small for the tightly bound 
nuclei up to ^^iVe and is getting larger with increasing number of neutrons. In this example the oscillator model would 
predict stable isotopes up to ^^A^e, while the calculation using the PW basis predicts stable isotopes up to ^°Ne. 

IV. SUMMARY AND CONCLUSION 

It has been the aim of this study to establish a tool, which allows a microscopic description of nuclei closed to 
the neutron drip and nuclear structures in the transition from isolated nuclei to homogeneous rnatter based on a 
realistic NN interaction. A model space technique, the unitary-modcl-opcrator approach (UMOA) jl9|, has been used 
to disentangle the low- momentum parts from the high- momentum parts from the CD Bonn interaction^. This leads 
to a universal effective interaction V/otufc, which is essentially independent on the underlying realistic interaction V , if 
the cutoff A is chosen around A = 2 fm~^. 

The Hartree-Fock (HF) calculations using Viowk do not give any saturation in nuclear matter and also HF calculation 
for finite nuclei lead to nuclear systems of high density and large binding energies. It has been suggested to compensate 
this feature by adding three-body forces [l3|. Instead we supplement the Viowk interaction by a density-dependent 
contact interaction and adjust the parameters in such a way that HF calculations for infinite matter reproduce the 
empirical saturation point for symmetric matter and the symmetry energy. 

This interaction model is than used to evaluate the properties of finite nuclei. For nuclei with closed shells for 
protons and neutrons we obtain a very good agreement with the empirical data. Special attention is paid to the 
technique for solving the HF equations for finite nuclei. The standard technique, expanding the single-particle wave 
functions in a basis of appropriate harmonic oscillator states is compared to a technique in which these states are 
expanded in a basis of plane waves. Rather close agreement between these approximation schemes is obtained for 
nuclei, which are tightly bound. The plane wave basis is favorable, however, for calculating loosely bound systems, 
since it allows for a better description of the asymptotic part of the nucleon wave functions. 

These tools, the interaction model and the technique for solving the HF equations, can now be used to describe 
the structure of loosely bound nuclei, like halo nuclei or nuclei at the edges of the valley of stability, or the nuclear 
structures in the so-called pasta phase of neutron stars. Another perspective is to consider a many-body approach 
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FIG. 7: (Color online) The neutron wave function of lpi/2 (dark lines) and the proton wave function of of ^"Ca (light 
lines) using two different basis systems, the plane wave basis (solid lines) and the harmonic oscillator basis (dashed lines). 
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FIG. 8: (Color online) The Fermi energies for the neutrons in the isotopes of Ne. Results of HF plus BCS calculations using 
the plane wave basis (PW) and the harmonic oscillator basis (HO) are shown for various even numbers of neutrons. 



for describing the nuclear systems, which goes beyond the HF approximation. 
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